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1 Introduction 

Understanding physics and stabilization of the Resistive Wall Mode (RWM) is 
an important task for the successful operation of the present day and the future 
fusion devices. Converted from the ideal MHD mode in the presence of the wall 
with finite conductivity, and growing on the time scale of the magnetic field 
penetration time through the wall, RWM sets pressure limit in the advanced 
scenario of the tokamak device [H HI |3] . It also limits the discharge duration 
of the Reversed Field Pinch (RFP) device HE |6]. Control of RWM growth is 
necessary for the discharge times longer than the wall time. 

Two general methods are known for the RWM stabilization: active control 
and stabilization via the mode resonance with continuum spectra or particle 
motions. For the active control the perturbed magnetic field measured by the set 
of sensors facing the plasma surface is used to generate the control signal by the 
set of active control coils. Such mechanism is extensively studied theoretically 
[3 IHl m Uni m HH and is successfully applied in the present day fusion devices, 
allowing operation with the plasma pressures above no- wall limit for tokamaks 
[TSmHHS] and resulting in substantial increase of the plasma discharge duration 
in RFPsdillT] . 

RWM suppression via the mode resonance with continuum spectra or par- 
ticle motions is another possible control mechanism (also called 'rotational sta- 
bilization'). It was experimentally observed for the tokamak configuration |18] 
that RWM is stabilized when the plasma rotation frequency is sufficiently high. 
Several mechanisms were proposed that could explain the stabilization. In ideal 
MHD description the dissipation caused by the RWM resonance with Alfven 
[TOl or sound HU [^^ continuum spectra could lead to the RWM stabiliza- 
tion. The dissipation strength depends on the plasma rotation frequency. The 
plasma rotation frequency n of the order of sound frequency is needed for 
the RWM stabilization by the resonance with continuum spectrum. Another 
mechanism that could lead to the RWM suppression in the region ft ^ ojg is 
ion Landau damping. Accurate description of this mechanism requires kinetic 
treatment of the ion motion parallel to the magnetic field [23] • In ideal MHD 
description this mechanism is modelled by adding dissipation term [M]. More 
physically consistent approach is used in semi-kinetic model |25J where damping 
term is calculated using drift kinetic energy principle for the large aspect ration 
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approximation. Mode resonance with circulating and trapped ions are included 
in the model. Calculated damping term are then included in the MHD descrip- 
tion. In order to explain observed RWM stabilization for the very low plasma 
rotation frequency [55] (fi <C Ws) other model was proposed [37] where the main 
dissipation channel is due to the mode resonance with the precession drift mo- 
tion of the trapped particles. This model in fact predicts the RWM stabilization 
even without plasma rotation. The accurate prediction of the RWM stabiliza- 
tion for the low plasma rotation frequencies is an important issue for the future 
reactor experiment ITER, where the highest plasma rotation in steady-state 
scenario is predicted to be in the sound range [28] . 

Previous studies of the rotational stabilization of the RWM were performed 
mainly for the tokamak configuration. Available results for the RFP [551 130] 
show that the plasma rotation frequency in the range of the Alfven frequency 
(fi ~ Wa) is needed for the RWM stabilization by the mode resonance with con- 
tinuum spectrum. Such high plasma rotation is not observed in the present day 
RFP devices, therefore rotational stabilization mechanism was not considered 
as realistic for RFP. On the other hand mode resonance with particle motions 
occurs for the low plasma rotation frequency and could be considered as a pos- 
sible RWM stabilization mechanism also for RFP configuration. In this work 
numerical studies of the effect of kinetic resonances on the RWM stability are 
performed for reversed- field pinch (RFP) configuration. Particularly, parame- 
ters corresponding to the RFX device [3T] is used. In RFP the values of poloidal 
and toroidal components are of the same order of magnitude. Toroidal compo- 
nent of the equilibrium magnetic field changes sign at the plasma edge. Several 
Fourier harmonics with different poloidal and toroidal mode numbers are seen 
in the RFP mode spectrum. Resonant harmonics with poloidal mode numbers 
m — 0,1 are usually linearly stable. Several non-resonant harmonics with dif- 
ferent toroidal mode numbers (positive and negative according to the sign of 
the safety factor q of the resonant magnetic surface) are unstable growing on 
the times comparable with the equilibrium magnetic field penetration through 
the resistive wall and therefore are classified as RWMs. Experimental studies 
of RWM control in RFP have shown possibility of simultaneous suppression of 
several unstable harmonics with different toroidal mode numbers using active 
feedback [HHT]. 



In this section physical model is briefly described that is used in the present 
work. More complete derivation can be found in [35]. Stability of RWM for 
these studies is determined by solving numerically the system of MHD equations 
with toroidal plasma rotation; 



p{-iuj + inn)v = -V-p+jxB + JxQ - p[2nZ x v + (v • Vr2)i?^ V.^] (2) 



2 Model 



{-iuj + mn)^ = V + (C • Vn)R^V(f> 



(1) 



{-ioj + inn)Q V X (v X B) + (Q • Vn)R'^\/(l) 
(—ico + inn)p = —V ■ VP 
j = VxQ 



(3) 
(4) 
(5) 



2 



whore uj is the complex mode frcqiiency (lu = i"/ — uir where 7 - mode growth or 
damping rate and LOr - mode rotation frequency), B,J,P- equilibrium magnetic 
field, current density and pressure respectively, Z - unit vector in the vertical di- 
rection, p - plasma density, - plasma rotation frequency in toroidal direction (j), 
C)V,j,Q - plasma displacement, perturbed velocity,perturbed current, perturbed 
magnetic field respectively, p - pressure tensor. 

Kinetic terms are included into the MHD equations via the pressure tensor 
components. Pressure tensor is defined as 

p = Ip + P||bb + p_L(I-bb) (6) 

where p is the scalar fluid pressure perturbation, p|| , pj_ are the particle parallel 
and perpendicular components of the kinetic pressure perturbations, b = B/i?, 
i? =1 B I, I is the unit tensor. For the particular perturbation with the mode 
number n, the parallel and perpendicular components of the pressure tensor are 
written as: 

p^^e-'"t+in<P J J dedKMvjfl (7) 

e.i 

^^6--*+-*=^! JdedA^Mvlfl (8) 

e.i 

Here the summation is over the electron and ion plasma components, integration 
is carried out over the particle energy e and pitch angle A, M is the particle 
mass, f II . v±_ are the parallel and perpendicular velocity components respectively 
with respect to the equilibrium magnetic field, is the perturbed particle 
distribution function. It is derived as: 

fl = -/°efce-'"*+™^ ^m^mi Am;e-™^(*)+'™<'^>+''"''* (9) 

where subscripts n,m,l mark Fourier components along toroidal and polidal 
angles and along bounce particle orbit, - derivative of the equilibrium dis- 
tribution function (taken to be Maxwellian for thermal particles) with respect 
to the particle energy e, ej, - kinetic energy of the particles, Hmi - perturbed 
Lagrangian component in Fourier space, \„d mode-particle resonance operator, 
4>{t) = (l){t)— < 4> > t - the periodic component of the particle velocity in 
toroidal direction, < .. > means average over the particle bounce period. The 
mode-particle resonant operator is calculated as: 

^ ^ _ n{uj^N + (efc - i/2)uJ.^.T + - co 

nu>d + nfl + [a{m + nq) + l]LOb — ii'eff — ^ 

where uj,,j\j,uj^t - diamagnetic frequencies associated with the density and tem- 
perature gradients respectively, ojb - bounce frequency, uid - bounce averaged 
magnetic drift frequency, Veff - eff'ective collision frequency, a =0 for the 
trapped particles, a =1 for the passing particles. The energy transfer between 
the mode and the particles is described by the imaginary part of the resonant 
operator. It can be seen that in case of finite mode growth (damping) rate and 
collision frequency energy is transferred between the mode and the particles for 
all plasma rotation frequency values. Mode-particle interaction is enhanced in 
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the resonant regions, where the real part of the resonant operator is equal to 
zero. The collisionallity effect on the RWM stability (via collision frequency 
term ii'eff of the resonant operator) is not considered in this work. 

Using dni) and (fTUj) to calculate the pressure tensor components and sub- 
stituting ^ in ([2]), the self-consistent formulation is obtained, that allows to 
study RWM stability including the effect from kinetic resonances. 

System of ideal MHD equations where the ion Landau damping is modelled 
by the parallel viscosity can be obtained from the equations (P) - (O by using 
only scalar pressure part of the pressure tensor adding the viscous stress 
tensor V-II = K\k\\Vthi\p{'^ • bb) in eq. ((!]), where k is the coefficient determining 
the damping strength, fc|| = (n — ^)/R, vthi is the ion thermal velocity and 
adding the term — 5/3-PV • v in the eq. (0]). 

3 Numerical details 

System of equations (HJ - ([5]) is an eigenvalue problem that can be written in 
the matrix form as 

-foBX = AX (11) 

where 70 is the complex eigenvalue, A,B- MHD operators, X - eigenvector. 
In this work the eigenvalue problem ([TT|) is solved numerically using MHD sta- 
bility code MARS-K in toroidal coordinates (s,x,0), where s = \J~^~~^ 
is the normalized poloidal flux {ipQ is the poloidal flux value on the axis), x is 
the generalized poloidal angle and (j) is the geometrical toroidal angle. Equilib- 
rium quantities are obtained using equilibrium solver CHEASE [3.3 . Note that 
eigenvalue problem is non-linear as resonant operator ()10p includes eigenvalue 
70 and therefore A — ^(79). Outer iterative loop is required in order to obtain 
the solution of the eigenvalue problem. 

4 Results 

4.1 Equilibrium 

Equilibria closely modelling the RFP device RFX ^31j is used. Circular plasma 
cross-section is taken with major radius Ro — 2.0 m and aspect ratio e — 
0.23. The RFP equilibrium parameters are F = ||^=-0.06, 6 = ^§^ = 1.41, 
where Bg{a),B^(a) are the values of toroidal and poloidal components of the 
equilibrium magnetic field at the plasma edge and < . . > means averaging over 
the plasma column. The plasma current value is /p=1.6 MA, on-axis toroidal 
field i?o=1.53 T. The poloidal beta value is (3poi=i%, the on-axis values of 
electron and ion temperatures and densities are: Teo=l kEv, Tio=400 eV, UgQ = 
nio=2.0- 10^^ m^'^. The pressure profile p{r) is given by the expression p{r) — 
n{r) * T(r), where n{r) = no{l - {r/af), T{r) = To{l - [r/af). The on-axis 
and edge safety factor values are: go=0.161, qa=-0.01. In Fig. [T]the poloidal 
and toroidal components of the magnetic field, pressure profile, and safety factor 
profile are shown for the described equilibrium as a functions of the normalized 
poloidal flux. 
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Figure 1: Equilibrium quantities as a functions of the normalized poloidal flux, 
a) normalized equilibrium magnetic field components (solid -toroidal component, 
dashed - poloidal component), b) normalized pressure, c) safety factor 



4.2 RWM stability without rotation 

The RWM spectrum in RFP is usually characterized by the presence of several 
unstable RWMs with different toroidal mode numbers n. Fourier harmonics with 
different helicities (positive and negative toroidal mode numbers) are visible in 
the spectrum due to the edge reversal of the equilibrium toroidal field. In 
Fig. [2] the spectrum of unstable RWMs is shown calculated for the equilibrium 
described in the previous section . The poloidal structure (different poloidal 
Fourier harmonics) of the normal displacement ^„ of the most unstable toroidal 
harmonic (n=-6) are shown in Fig. [S] It is observed that the harmonic with 
m=l is dominant. This result points to the weak poloidal coupling in the 
RFP (differently from the tokamak configuration). In the further discussion 
the stability of the Fourier harmonic with n=-6 will be discussed if not stated 
otherwise. 

4.3 Ideal MHD stability with plasma rotation. Contin- 
uum resonances 

Plasma rotation opens the possibility for the interaction of the static non- 
resonant RWM with the stable waves in the plasma. In ideal MHD such in- 
teraction is due to the mode resonance with continuum spectrum. Two con- 
tinuum spectra are known to be resonant with RWM: Alfven continuum and 
sound continuum. The condition for the resonance appearance can be written 
in general as nD, + w,. = lOc where Q, is the plasma rotation frequency, uj^ is 
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Figure 2: Spectrum of unstable RWMs. Black - n < 0, grey - n > 0. 
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Figure 3: Poloidal harmonics (m=-5..5) of the RWM with n~-6. Dominant 
m^l poloidal harmonic is shown in magenta solid line. All other harmonics are 
not visible. 
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Figure 4: Mode resonance with the Alfven continuum. Complex eigenvalue 
(mode growth rate - and mode frequency uJrTw) dependence on the plasma 
rotation ^l/uaP normalized to the on-axis value of Alfven frequency. 

the RWM rotation frequency and Wc is the continuum frequency. For the static 
RWM (ojr = 0) this condition becomes nVt = lOc- RWM resonance with the 
Alfven continuum spectrum appears when the condition \nQ,\ = \ujca\ = |fc||^^a| 
is satisfied. Here Va = B /{nopY^'^ is the Alfven velocity and fc|| = {m/q ~ n) is 
the parallel component of the RWM wave vector. In Fig. 2] the mode growth 
rate 7 and mode rotation frequency ojr (both normalized by the wall time r^) 
dependence is shown on the plasma rotation frequency f2 (normalized by the 
Alfven frequency on the plasma axis wj^). The pressureless equilibrium with the 
equilibrium parameters F, 6 given in sec 14. H is used. 

Because of the non-resonance nature of the RWM in RFP, no resonance be- 
tween the mode and the Alfven waves is possible at vanishing or slow plasma 
rotation. The obtained value of the plasma rotation frequency that starts to 
give the resonance with Alfven wave (seen at peak of the growth rate •yT^) is 
n/LU° w0.06. This critical value agrees very well with the analytical prediction 
that is calculated using the resonance condition given above. The RWM growth 
rate "/Tw starts to decrease for plasma rotation frequencies 57 > LUca accompa- 
nied by the finite mode frequency appearance. This points to the damping effect 
appearance from the mode-continuum resonance, although no compete RWM 
stabilization is seen for the studied O range. The obtained value of the plasma 
rotation frequency for the mode resonance with the Alfven continuum is much 
less that the one obtained in the previous studies [30] due to the fact that the 
internal RWM with n = — 6 is used in the present studies with the resonant sur- 
face closer to the plasma (and therefore smaller A;||). It should be noted that the 
experimentally observed plasma rotation frequency in RFX is still much smaller 
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Figure 5: Complex eigenvalue (mode growth rate - and mode frequency 

UJrTw) dependence on the plasma rotation il/wa normalized to the Alfven fre- 
quency. Blue (and x )- k—0.0, green (and a) - k—0.5, red (and D) - k — 1.0, 
cyan (and ()) - k~ 1.5 

than the value obtained here for the mode resonance with Alfven continuum (fi 
value is of order of lQ~^Wa [31] )■ 

The resonance of the static RWM with sound continuum appears when = 
Wcs = k\\Vs where Vg = {p/ pY^"^ is the sound velocity. The behaviour of the 
complex eigenvalue in the plasma rotation frequency range corresponding to the 
resonance condition is shown in Fig. [5] (line with x) where an equilibrium with 
finite pressure is chosen. The obtained plasma frequency value that corresponds 
to the resonance is ^/oj^ ~ 0.004 and also agrees well with the value given by 
the analytical predictions. Note also that lOcs ^ca- 

According to the studies performed in |23j the resonant behaviour seen at Ucs 
is not physical due to the fact that the ideal MHD model is unable to describe 
accurately the particle motion along the field lines. A reasonable damping model 
can be introduced [24j to remove the unphysical resonant behaviour observed 
in Fig. [S] by damping the parallel sound wave in the ideal MHD model. This 
parallel wave is damped physically by ion Landau damping mechanisms. In this 
work viscous damping model is used (see Sec. [2|) where a numerical coefficient 
K, is introduced to measure the strength of the parallel sound wave damping. 
In Fig. [5] the complex eigenvalue behaviour is shown for the different values 
of the sound wave damping coefficient. It is seen that the the large k values 
(strong damping) removes the resonant behaviour of the RWM growth rate. At 
the same time viscous damping results in RWM suppression, with the larger 
suppression rate for the larger k values. 



8 



4.4 Kinetic resonances 



A more physically consistent description of RWM stability for the low values of 
the plasma frequency is obtained considering the mode interaction with plasma 
particle drift motions. Mode particle interaction is introduced in the MHD 
model through the pressure tensor ^ and is included into the resonant operator 

^ _ 7l(w*Ar + (cfe - 3/2)w*t + fl) - CJ 

^ ^ ^ nuJd + nVl + \a{m + no) + — iVeff — ^ 

where sum is over the poloidal Fourier harmonics m, bounce harmonics I, par- 
ticle fraction a, and particle species (e, i). It is seen that resonance with several 
particle motions is possible depending on the value of the plasma rotation fre- 
quency. In these studies the mode resonance with precession drift and bounce 
motions is considered. 

Full picture of the mode-particle resonance is complex due to the fact that 
contribution from the both trapped and passing fractions of the different plasma 
species (electrons and ions) should be taken into account. Moreover, particular 
particle motion (characterized by the motion's frequency value) have in general 
complex dependence on the pitch angle A and particle kinetic energy on each 
flux surface. Total effect of the RWM-particle resonances on the mode stability 
is obtained numerically considering mentioned dependencies. Kinetic frequen- 
cies averaged both over e and A are discussed in this section. Such simplified 
approach allows to make some qualitative predictions for the mode-particle res- 
onances. The average kinetic frequencies for the ion plasma component are 
shown in Fig. [B] as functions of the normalized poloidal flux. The frequency 
for the passing particles ujp is equal to Wb when a = 1 in (jlOp . Note that LUp 
and uJb are scaled (by factors 0.01 and 0.1 respectively) to be able to compare 
profiles of different kinetic frequencies on one figure. The value of is the sum 
of the diamagnetic frequencies due to the density and temperature gradients. 
The plasma rotation frequency is not shown on the figure above, but comparing 
the experimental value (fi ~ lO^'^Ua 34 ) with the calculated results it is seen 
that fl ujd- The following approximate frequency ordering is obtained for the 
studied case: 

~ Wrf w* <C Wf, ~ ojca ~ ojf < Wp < a;° (12) 

where cof ~ 0.1w° is the central ion sound frequency. It is seen that the ion 
bounce frequency value is in the range both with the sound frequency and also 
with the frequency of the mode resonance with Alfven continuum. 

In the present work the plasma frequency range ft < ujf is studied. Two 
specific frequency subregions can be deduced following the frequency ordering 
obtained above: D, ^ ujd and uj ^ ujf,. In each of this frequency region the plasma 
rotation frequency is close to the certain kinetic frequency (precession drift of 
bounce frequencies) and therefore the mode-particle resonance is expected to 
affect the RWM stability. It is useful for the further studies to have qualitative 
estimation of particle fractions and particle species that possibly could be in- 
volved in the mode-particle resonance in the mentioned frequency regions. The 
following qualitative analysis aims the estimation of the dominant terms of the 
resonant operator A for the two frequency subregions for different values of a 
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Figure 6: Ion kinetic frequencies averaged over the pitch angle A and particle 
kinetic energy e^. Solid line -ujb*0.1, dotted line - ujp*0.01, dashed line -Ud, 
dashed-dotted line - w*. All frequencies are normalized to the central Alfven 
frequency wafaew^. 

(particle trapping) and I (bounce harmonic number). In the subregion i7 ^ 
the following conditions could lead to the mode-particle resonance: 

- a = 0. In this case A ~ — , nltw > — • For I = mode resonance with 
precession drift motion has dominant effect. Both electrons and ions contribute. 
For / ^ bounce motion could contribute to the mode-particle resonance for the 
slow rotation when ujd ~ {l/n)uJb- As ujd/'^b ~ 0.02 across most of the plasma 
column (see Fig. [B]) condition I <C n should be satisfied that is not the case in 
the present studies {min\{l/n)\ ~ 1/6). Note that ood ^ is seen locally in the 
plasma core or the plasma edge regions and therefore bounce motion can affect 
the RWM stability in these regions. 

- a = 1. In this case resonant operator has general form ([TU]). Passing 
particles can contribute to the mode-particle resonance for the slow plasma 
rotation when {m + nq -\- l)/n ~ ujd/oop. As \uJd/ojp\ < 10~^ in the present 
studies the contribution from the passing particles can be expected when m = 
— /, l^l 1.0 that is satisfied locally near the resonant surface. 

The main effect on the RWM stability in the low rotation frequency sub- 
region is expected to come from the mode resonance with the precession drift 
frequencies. The effect from the bounce motion or passing particles is local. 

Similar analysis for O ~ gives: 

- a = As from the frequency ordering iOd <S oJb mode resonance only with 
bounce motion is possible, i.e. when I ^ 0. Then A ~ (;/n)X^o • mode- 
particle resonance condition is satisfied when {l/n)ujb ~ —57. For present studies 
(n — —6) bounce harmonic with I — 6 should have the dominant contribution. 
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Note also that this condition could give the contribution to the RWM stability 
in the intermediate values of the plasma rotation frequency (i.e LOd ^ ^ < Wfe)- 
This happens when l/n < 1. Only ion component is involved as the electron 
bounce frequency is much higher. 

- a = 1 In this case A ~ 77 — , "i!7tP 1 — to- Passing particles can contribute 
to the mode resonance with the bounce motion when (m + nq + l)/n ^ ujb/ujp. 
As jwb/wpl ^ 10^^ in this studies the contribution of the passing particles can 
be expected for I = — m, |g| < 0.1 that points to the possibility of the global 
effect from the passing particles (note that \q\ < 0.16 in the present studies). 

Both trapped and passing particles could contribute to the mode-particle 
resonance for the plasma rotation in the order of ion sound frequency (note 
that ujs ^ ujb from the frequency ordering (jl2p ). Trapped particle contribution 
comes from the bounce motion, passing particles effect could be global (note 
that this is different from the case of the low rotation frequencies $7 ^ w^)- 

4.5 Comparison of the viscous and kinetic damping mod- 
els 

It was mentioned above that the reasonable damping model can be used in ideal 
MHD description to model physical ion Landau damping. Viscous damping 
model is used here with the damping coefficient k=1.0 (strong damping) for the 
comparison with the MHD description including kinetic effects. The results of 
comparison are shown in Fig. [71 It is seen that kinetic damping model gives 
slightly lower growth rate for the low plasma rotation frequencies (fl ^ uJd) while 
the higher suppression rate is obtained for D, cob using the viscous damping 
model. Increase of the RWM growth rate for > 5 • 10~^cij° is explained by 
the resonant RWM behaviour in vicinity of the mode resonance with Alfven 
continuum (see frequency ordering (I12[) '). 

4.6 Different kinetic resonances 

Several kinetic frequencies are included in the resonant operator (|10l) and there- 
fore can contribute to the total effect on the RWM stability. In the previous 
section 'full' model is used taking into account mode resonance with precession 
drift and bounce motions for trapped and passing particles. In order to evaluate 
the importance of the particular kinetic resonance and particle fraction on the 
total effect the following studies are performed. It can be seen from the def- 
inition of the resonant operator that the frequencies of the kinetic resonances 
are included in additive manner, giving the possibility to exclude all but one 
resonance from the calculations. Particle fractions (trapped or passing) can be 
also treated separately. In Fig. |S]the 'full' model (blue curve) is compared to 
the 'only bounce' model (red curve) when only bounce frequencies are taken 
into account and with 'only precession drift' model (green curve) when only 
precession drift frequencies are taken into account. Both trapped and passing 
particles are included. Two plasma rotation frequency regions can be distin- 
guished separated by the plasma rotation frequency value ^/uj^ ~ 10~^. For 
the region fi/wj^ < 10~^ eigenvalue calculated for the 'full' model case is closely 
follows that calculated for the 'only precession drift' model. Note that the mode 
growth rate for these two cases is slightly less than that for the 'only bounce' 
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Figure 7: Complex eigenvalue (mode growth rate - "fT^ and mode frequency 
UJrTw) dependence on the plasma rotation fi/wa normalized to the Alfven fre- 
quency. Solid line with x- fluid damping (sound wave damping with k — l.Q), 
dashed line with □ - model with kinetic resonances included. 

model. For the region fi/a;" > 10"'^ situation becomes opposite i.e. the eigen- 
value trace corresponding to the 'full' model follows that for the 'only bounce' 
model. Resonant behaviour of the RWM growth rate is seen for the 'only preces- 
sion drift' model in this frequency range. It is attributed to the mode resonance 
with the Alfven continuum spectrum. 

In Fig. ini the same resonances are compared where only trapped particles 
are included. It is seen that the mode growth rate behaviour is different in the 
second plasma rotation frequency region Vl/uj'^ > 10^^ for the 'full' and 'only 
bounce' models. Now similar resonant behaviour is seen for all three models at 
ri/a;" w 0.05. The growth rate for the 'full' model is almost identical with that 
for the 'only precession drift' model for the whole range of the studied plasma 
rotation frequencies. The effect of the mode resonance with different particle 
fractions for two frequency regions (r^/w^ < 10^^ and JI/wq > 10~^ ) can be 
seen. In the region fi/o;" < 10^^ the effect of the trapped and passing particles 
on the RWM stability is observed. RWM is suppressed due to the resonance 
with the precession drift motion (x curve in Fig. [5]) and destabilized due to 
the resonance with the passing particles(o curve in Fig. [5]). This two effects 
almost cancel each other as it is seen from the behaviour of the 'full' model 
eigenvalue (□ curve in Fig. Note that there is no effect from the bounce 
motion in this frequency region (compare a curves in Figs |8] and |9l). In the 
region il/uj'^ > 10~^ only passing particles have the effect on the RWM stability. 
Indeed, similar eigenvalue behaviour for the three studied cases seen in Fig. [S] 
points (somewhat surprisingly) that the mode resonance with bounce motion 
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Figure 8: Complex eigenvalue (mode growth rate - -yTy, and mode frequency 
uirTw) dependence on the plasma rotation il/uja normalized to the Alfven fre- 
quency. Solid line with □ -full model, dashed line with x - resonance only with 
LOd is considered, dashed- dotted line with a - resonance only with cob is considered. 
Both trapped and passing particles are included 
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Figure 9: Complex eigenvalue (mode growth rate - -yTy, and mode frequency 
uirTw) dependence on the plasma rotation il/uja normalized to the Alfven fre- 
quency. Solid line with □ -full model, dashed line with x - resonance only with 
LOd is considered, dashed- dotted line with o - resonance only with cob is considered. 
Only trapped particles are included 
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does not affects the RWM stabifity. Tfie resonant beliaviour and subsequent 
mode suppression seen for Jl/w^ > 5 • 10~^ is attributed to tlie mode resonance 
witli the Alfven continuum. On the other hand substantial mode suppression 
is observed in Fig. [5] (o and □ curves) in this frequency range that effectively 
prevents the mode growth due to the resonance with Alfven continuum. Note 
that although the frequency of the passing particles is in general much higher 
than the frequencies of the bounce and precession drift motions, the possibility 
of the mode resonance with the passing particles for low frequencies exists and 
the conditions for such interaction are qualitatively estimated in Sec. 14.41 (cases 
a = 1 for bounce and precession drift frequency ranges). 

4.7 Equilibrium plasma pressure 

In the previous studies of the kinetic effects for tokamak configuration [35j it 
was shown that equilibrium pressure affects the mode interaction with the par- 
ticle motions. Both plasma pressure profile and absolute value can change the 
behaviour of the complex eigenvalue. In order to investigate the effect of pres- 
sure on the mode resonance with the plasma particles for RFP configuration 
the studies of the plasma pressure profile pressure absolute value (characterized 
by the f3p parameter) are performed below. 

4.7.1 Pressure profile 

Pressure profile shape affects the value of the pressure derivative and therefore 
changes the values of frequencies and also (according to [^) affects the 
precession drift frequency term of the resonance operator (jlO|) . Two pressure 
profiles are compared in the present studies: the one used for the studies above 
(see Sec l4.ip and more flat one that is characterized by the following density and 
temperature profiles: n'^%r/a) = n^^^il - {r /a)^2),T^'''(r/a) = T^^l - (r/a)^). 
The results of comparison are shown in Fig. [TUl It is see that by changing 
pressure profile the diamagnetic frequency profile is changed (dashed-dotted 
line in Fig. [TObV Note also the slight change of uJd (dashed line) and Wp (solid 
line) profiles. The RWM is suppressed more (Fig. [TUb ) in the region J7 ~ for 
the case of the fiat pressure profile, but for ^ cut more peaked profile gives 
better suppression. RWM destabilization is seen for > 5-10~^a;° in both cases 
attributed to the effect from the mode resonance with the Alfven continuum. 



4.7.2 EflFect of /3p 

The absolute value of the plasma pressure (characterized usually by the beta 
value) is one of the most important parameters in the fusion research. Operation 
with the high beta value is necessary for the efficient operation of the fusion 
reactor. In this work poloidal beta /?„ = ^^^jf^Vtot is used for the pressure 
value characterization, where < p > is the equilibrium plasma pressure averaged 
over the plasma volume Vtot, Ip - plasma current. It is seen from the frequency 
ordering (jl2p that for n = — 6 Fourier harmonic the effects from the kinetic and 
continuum resonances occurs for the same plasma rotation value lo ^ ujf. In 
order to separate these effects another unstable Fourier harmonic is chosen with 
n=-5 for which the mode resonance with the Alfven continuum is shifted towards 
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Figure 10: Results of the pressure profile studies, a) normalized pressure radial 
profiles; h) kinetic frequencies (solid - Wp * 0.01, dashed - lOd, dashed-dotted - 
CO*), cob profiles are not shown as no effect on RWM stability is found from the 
mode resonance with the bounce motion ; c) RMW growth rate dependence on 
the plasma rotation frequency; d) RWM rotaMon frequency dependence on the 
plasma rotation frequency. Black - n{r) — no(l — {r / of' ,T{r) = To(l — {r/a)^; 
grey - n(r) = no(l - {r/ay^,T{r) = To{l - {r/af; 
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Figure 11: (Colour online) Results of the pressure value studies for the Fourier 
component with n = —5. a) normalized pressure radial profiles; b) kinetic fre- 
quencies (solid - cjp * 0.01, dashed - LOd, dashed-dotted - uj^,); c) RMW growth 
rate dependence on the plasma rotation frequency; d) RWM rotation frequency 
dependence on the plasma rotation frequency. Blue - Pp=5%; red - j3p=15%, 
green -/3p=25%; 

higher plasma rotation frequencies. Indeed using analytical expression (see Sec. 
14. 3p the mode resonance with Alfven continuum appears for this harmonic, when 
n PS 0.2a;°. The results of the studies for three /3p values (/3j,=5%,15%,25%) are 
shown in Fig. [TT] 



It is observed (Fig. [TTk ) that the values of kinetic frequencies depend on poloidal 
beta. Dependence for different kinetic frequencies can be approximated as: 
Wp(/3p) ^ a/^, ^jJdiPp) ^ Pp, 0J^,{f3p) ^ (3p that is similar as for the large aspect 
ration scaling for tokamak equilibrium [35] • Change of the kinetic frequen- 
cies magnitude causes also the change of the plasma rotation frequency of the 
mode-particle resonances. The consequence of such resonance shift is that the 
stabilizing effect from the kinetic resonances in the range f2 ^ cjf, is seen only 
for the low value of poloidal beta (blue and red curves In Fig. ITTb ) for the 
studied range of plasma rotation frequencies. The mode destabilization seen for 
the > 0.01 is similar for all cases and can be attributed to the effect of the 

mode resonance with Alfven continuum. Note that no such substantial shift of 
the kinetic frequencies is seen for tokamaks, due to the fact that the beta range 
where unstable RWM is observed is of the order of one percent that is much 
smaller than in RFP. 
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4.8 Mode resonance with electron and ion plasma compo- 
nents for Q ~ 

It was shown in [32j for tokamak that temperature difference between electron 
and ion components can sufficiently affect the RWM stability for the low plasma 
rotation frequencies (il ^ ujd) pointing to the different contribution from the 
mode resonance with electrons and ions. In general, contribution both from ion 
and electron components to the total effect is present at each flux surface due 
to the specific dependence of the precession drift frequency on the pitch angle 
for each fiux surface. In particular, ujd changes sign as a function of A (see Fig. 
1 in ^32) ) and therefore mode resonance both with electrons and ions is possible 
on the same fiux surface as ujd depends on the sign of the particle charge (note 
that such u)d{^) dependence makes possible RWM suppression with zero plasma 
rotation). However for the present studies it is seen (Fig. ^ that the average 
ion precession drift frequency is positive for all radial points. This means that 
precession drift frequency for particular plasma component has predominantly 
one sign as a function of A on each fiux surface and therefore the contributions 
from electrons and ions could be separated by the sign of the plasma rotation 
(i.e. mode resonance only with ion or electron component will determine the 
total effect on the RWM stability for the particular plasma rotation direction). 
In Fig. [12] the complex eigenvalue behaviour is shown as a function of the 
plasma rotation frequency in the low plasma rotation region (f2 iOd). The 
complex eigenvalue is calculated for the Fourier harmonic n — ~6 taking into 
account only the mode resonance with the precession drift motion. Two values 

of the ion- electron temperatures ratio Cp = ( C„ = 0.33 and C„ — 0.5) 

are compared for the positive and negative plasma rotation directions. First Cp 
value corresponds to the ion-electron temperature ratio used in the present work 
(Teo=l keV, Tio=0.5 keV), second Cp value corresponds to the case of equal 
central temperatures of electrons and ions. In this case = —uj^,ujd = —<jJd- 
Absolute value of the mode rotation frequency (|aJr|Tu,) is plotted in order to 
make comparison easier (mode frequency rotation changes sign when the plasma 
rotation direction is reversed). It is seen that the mode growth rate behaves 
differently for the positive and negative plasma rotation directions in the case 
when Cp = 0.33. For negative O (dashed line with 0) mode stabilization effect 
is observed for the lower values than for the positive rotation direction (solid 
line with □). This could mean that |a;^| < |aj^| and the mode resonance with 
electron component occurs for the lower J7 values. In the case when Cp = 0.5 
the same effect is observed for the positive and negative fl directions (solid line 
with <l and dashed line with > are overlapped). This is expected result as the 
kinetic frequency values are equal for Cp = 0.5 and therefore the effect from the 
different plasma components is independent on the plasma rotation direction. 



5 Discussion and conclusions 

The results of the present investigations show that the kinetic effects do not 
change the RWM stability substantially for the studied RFP equilibrium. An 
attempt to explain such result could be made by considering the dispersion re- 
lation for the RWM written in terms of the fiuid and kinetic energy components 
(see eq. 8 in [37]). This dispersion relation is used for the studies of the ki- 
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Figure 12: Complex eigenvalue dependence on the plasma rotation frequency for 
the Fourier component with n = —6. Solid line with □ - O > 0, Cj, = 0.33, solid 
line with < - Q > 0,Cp — 0.5, dashed line with <) - < 0,Cp = 0.33, dashed 
line with > - Q < 0, Cp = 0.5. Mode growth rate 7 for each curve is normalized 
to the value for O = 10~^. 
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netic effects via the perturbative approach. Quahtatively, the stabihty of RWM 
is determined by the relation between fluid and kinetic energy components. It 
was shown [36j that for the pressure driven RWM seen in tokamak configuration 
kinetic and fluid terms are comparable. Therefore changes in kinetic term can 
substantially affect the overall stability of RWM. The calculations made for the 
studied here RFP equilibrium show that for the current driven RWM fluid term 
is much larger that the kinetic term. Therefore qualitatively kinetic term does 
not play an important role for the RWM stability in RFP. More work should 
be done in order to make the general conclusion about the role of the kinetic 
energy on the RWM stability in RFP. 

It is also seen during the present studies, that the observed effect on the 
RWM stability from the kinetic resonances comes mainly from the mode reso- 
nance with the passing particles for the plasma rotation of the order of ion sound 
frequency. Surprisingly no effect from the trapped particles (bounce motion) is 
seen in this region. It is observed that also for the low rotation (of the order of 
precession frequency) mode resonance with the passing particles have quantita- 
tively similar effect on RWM stability as compared to the mode resonance with 
the precession drift motion. 

In conclusion the studies of the RWM interaction with particle motions was 
investigated for RFP configuration. The equilibrium close to the experimental 
was used. The effect from the kinetic resonances is small, but the eigenvalue 
behaviour for the low plasma rotation frequency is explained better by the 
model with kinetic resonances than the ideal MHD model. Two regions of 
stabilization was found similarly to the previous work[] that correspond to the 
different plasma particles. 
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